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ABSTRACT 

We critique the method of constructing extinction curves of lensing galaxies 
using multiply imaged QSOs. If one of the two QSO images is lightly reddened 
or if the dust along both sightlines has the same properties then the method 
works well and produces an extinction curve for the lensing galaxy. These cases 
are likely rare and hard to confirm. However, if the dust along each sightline 
has different properties then the resulting curve is no longer a measurement of 
extinction. Instead, it is a measurement of the difference between two extinction 
curves. This "lens difference curve" does contain information about the dust 
properties, but extracting a meaningful extinction curve is not possible without 
additional, currently unknown information. As a quantitative example, we show 
that the combination of two Cardelli, Clayton, & Mathis (CCM) type extinction 
curves having different values of Ry will produce a CCM extinction curve with 
a value of Ry which is dependent on the individual Ry values and the ratio of 
V band extinctions. The resulting lens difference curve is not an average of the 
dust along the two sightlines. We find that lens difference curves with any value 
of Ry, even negative values, can be produced by a combination of two reddened 
sightlines with different CCM extinction curves with Ry values consistent with 
Milky Way dust (2.1 < Ry < 5.6). This may explain extreme values of Ry 
inferred by this method in previous studies. But lens difference curves with more 
normal values of Ry are just as likely to be composed of two dust extinction 
curves with Ry values different than that of the lens difference curve. While 
it is not possible to determine the individual extinction curves making up a 
lens difference curve, there is information about a galaxy's dust contained in 
the lens difference curves. If the lens difference curve can be fit with the CCM 
relationship (regardless of the fitted Ry value), this implies that the dust along 
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the two slightlines can be described by CCM. In addition, the presence of the 
2175 A feature in the lens difference curve means that this feature is present 
along at least one of the two lensed sightlines. 

Subject headings: ISM:dust, extinction — galaxies:ISM — gravitational lensing 

1. Introduction 

The well-known Milky Way extinction relation of Cardelli, Clayton, & Mathis (1989, 
CCM) applies to a wide variety of Galactic interstellar environments. CCM extinction 
curves are described by one parameter, Ry, the ratio of total-to-selective extinction in the 
V band. The Ry value is a rough measure of average dust grain size, and therefore provides 
a physical basis for the variations in extinction curves. For the Galaxy, the often-quoted 
average Ry ~ 3.1, merely reflects the typical value for dust along diffuse sightlines. But 
the diverse environments within the Milky Way show a range of values, 2.1 < Ry < 5.6 
(Valencic, Clayton, & Gordon 2004). A Milky Way UV extinction curve is characterized 
by a "bump" at 2175 A and a rise in the far-UV. While CCM applies to most sightlines 
in the Milky Way and a few in the Magellanic Clouds, it does not typically apply in other 
extragalactic environments where extinction curves have been measured (Gordon et al. 2003; 
Gordon, Calzetti, & Witt 1997; Pitman, Clayton, & Gordon 2000). 

Traditionally, the pair method is used to determine extinction by comparing the spectral 
energy distributions (SEDs) from reddened and unreddened stars of the same spectral type 
(Massa, Savage, & Fitzpatrick 1983). The pair method has only been applied within Local 
Group galaxies such as the Milky Way, SMC, LMC, and M31 because it is necessary to 
use individual stars as point sources. Determining extinction in galaxies beyond the Local 
Group is generally a more complex problem. Any surface photometry of a galaxy contains 
light from a mixture of gas, dust, and stars, making it difficult to analyze the spectral energy 
distributions from these complex systems directly. To derive the SED of a galaxy requires 
a radiative transfer model which accounts for the physical properties of the dust, the stellar 
population, and the geometrical distribution of the dust within the galaxy (e.g., Gordon et 
al. 2001; Misselt et al. 2001). Applied to starburst galaxies, these models imply that these 
galaxies have dust with properties more similar to that in the SMC than to the Galaxy, 
having no measurable 2175 A bump (Calzetti, Kinney, & Storchi-Bergmann 1994; Gordon 
et al. 1997). The existence of the CCM relation in the Galaxy, valid over a large wavelength 
interval, suggests that the environmental processes that modify the grains are efficient and 
affect all grains. In other galaxies, environmental processing by radiation and shocks near 
regions of star formation or variations in galactic metallicity may be responsible for the large 
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observed variations in extinction curve properties from CCM (e.g., O'Donnell & Mathis 1997; 
Gordon & Clayton 1998). Certainly, the application of the CCM extinction law for Galactic 
dust is not appropriate for dust in starburst galaxies and it is unkown if it is appropriate for 
dust in more quiescent galaxies. 

In order to study extinction in galaxies beyond the Local Group, a new method has 
been developed. It uses multiply imaged QSOs produced by gravitational lensing to study 
the extinction properties of the lensing galaxy (Nadeau et al. 1991). The multiple images 
originate from the same QSO and thus have the same intrinsic SED. Therefore, this seems 
ideal for constructing extinction curves by the pair method. The different images are pro- 
duced by light traveling along different sightlines through the lensing galaxy, so the amount 
of reddening varies from image to image. The SEDs of the two images are then compared 
using the pair method to determine an extinction curve for the lensing galaxy. In this pa- 
per, we investigate the usefulness of this method in producing extinction curves of lensing 
galaxies. 



2. Background 

The literature currently contains ~25 multiply imaged QSOs, listed in Table 1, for 
which extinction curves have been calculated, using either spectroscopic or photometric data. 
Typically, a CCM relation was fit to these curves and a value of Ry found. Gravitational 
lensing is an achromatic process, but often the various images of a QSO are found to have 
different colors. The color differences are likely due to reddening by dust in the lensing galaxy 
(Yee 1988). Each image produced by a gravitational lens is also magnified by a different 
amount depending on the geometry of the source, and the lens, and possibly amplified due 
to microlensing effects. Early work by Nadeau et al. (1991) on the lens system Q2237+0305 
showed that each of four images from a single QSO source had different colors, and a method 
for determining the extinction was suggested. In order to find A(\)/A(V), the ratio of 
total extinction at a given wavelength, A, to the total extinction in the V-band, directly 
from the ratio of fluxes from two images, it is necessary to know the amplification ratio. 
Since the amplification ratios are uncertain, they assumed a value of A(X)/A(V) at the 
central wavelength of the K-band. Because the images only differed in magnitude by a 
small amount, Nadeau et al. assumed that the extinction at any given wavelength differs 
only slightly from the average Galactic extinction law. Nadeau et al. (1991) determined 
an extinction law for the lensing galaxy using the assumption that A(K)/A(V) = 0.11, 
the Galactic value (Clayton & Mathis 1988). The resulting extinction curve is remarkably 
similar to the Ry = 3.1 Galactic extinction curve, for 0.5 < A -1 < 2 /im" 1 . 
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Falco et al. (1999) built upon the work by Nadeau et al. (1991) and completed the first 
survey of extinction properties of 23 lensing galaxies. Using photometric data for multiply 
imaged QSOs, Falco et al. determined the differential extinction, AE(B — V), for each set of 
images. The process of determining differential extinction depends on several assumptions. 
First, they assumed that the source spectrum (from the QSO) is identical for each image. 
Second, they assumed that the spectrum does not show significant time variability over a 
span equal to the time delay for the images. Finally, they assumed that the magnification, 
AM, by the lens does not depend on wavelength or time. In order to compare the images, 
the magnification of one image was set to one, and the extinction of the bluest image set to 
zero. In several cases, Falco et al. use radio data to determine AM. They use some spectro- 
scopic redshifts, but also use a technique suggested by Jean & Surdej (1998) to determine 
photometric redshifts for several lenses. It should be noted that microlensing effects are not 
considered here, but Falco et al. suggested that these effects are small compared to the 
mean magnification, which is constant and achromatic. In many cases for which radio data 
were unavailable, the differential magnification, AM, was treated as a free parameter when 
determining the best fit extinction curve. Extinction curves have also been calculated from 
multiply imaged QSOs by Toft, Hjorth, & Burud (2000), Motta et al. (2002), and Munoz et 
al. (2004). 

Motta et al. (2002) use spectroscopic data to determine differential extinction in the 
lens galaxy SBS 0909+532. They suggested that spectroscopy is more suitable for the study 
of extinction in lenses because broadband photometry sometimes makes extinction and grav- 
itational microlensing indistinguishable. Microlensing causes the flux for a single image to 
increase inhomogeneously across the wavelength range of an individual filter, affecting the 
color of the individual image. Using spectroscopy allowed them to overcome many of the 
challenges presented by photometric methods. Wucknitz et al. (2003) used spectra of the 
doubly imaged QSO HE 0512-3329 to distinguish between the effects of extinction and mi- 
crolensing. Because extinction affects the emission lines, but microlensing does not, the 
emission lines can be used to separate these effects. Then, the continuum, which is affected 
by both microlensing and extinction, can be corrected for extinction and the microlensing 
effects may be studied directly. While extinction curves were not the focus of this paper, 
they fit the CCM relation to their data. They treated Ry as a free parameter, and fit their 
data using several values of AE(B — V) and AA(V). They suggested that the measured 
extinction curve could be the result of dust with two different values of Ry along the two 
sightlines. 

Each of the papers cited above uses roughly the same method for obtaining an extinction 
curve by the pair method. When calculating the extinction curve, they consider Ry to 
be a free parameter as well as the magnification, if it is not known. For several of the 
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lensing systems, observations in the radio, where extinction is negligible, allow for a direct 
measurement of the differential magnification (AM). The magnitudes of two images, labeled 
A and B, at wavelength A, are compared using the equation 

m B (\) - m A (\) = AM + AE(\)R(-^) (1) 

1 + Zi 

where AE is the differential reddening, and R is the ratio of total-to-selective extinction for 
the given wavelength and lensing galaxy redshift (Motta et al. 2002). This is similar to the 
traditional pair method where m B (\) — m^(A) = AE(X)R(X) = A(X). The lensing galaxy 
extinction curves are typically presented unnormalized making it difficult to compare these 
results to other extinction studies. Extinction curves should to be normalized to E(B — V) or 
A(V) so that only the wavelength dependence, not the amount of extinction, is compared. In 
this paper, all plots are presented in rest wavelength space and with the extinction normalized 
to A{V). It is straightforward to determine the normalized extinction using the following 
equation 

MX) = m B {\) - m A {\) - AM 

A(V) AE(B-V)R V [ ' 



3. Discussion 

Figure 1 shows a normalized extinction curve for one of the lenses, B 1600+434 (Falco 
et al. 1999). Using the photometric data in Falco et al., as well as their values for Ry, 
AE(B-V), and redshift, we plot the data and the fit which has an extremely low R v = 0.92. 
As pointed out by Falco et al. and shown in the figure, the small number of photometric 
points covering a limited wavelength range can be fit with CCM curves showing a very large 
range of Ry values. 

The results for the extinction curve given in Motta et al. (2002) for SBS 0909+532 differ 
greatly from those given by Falco et al. (1999) for the same QSO. In fact, the differential 
extinction is similar (both reported values around AE(B — V) = 0.2) but the inferred values 
for Ry are different. Falco et al. (1999) used an estimated redshift of z — 0.60 for the 
lensing system and get Ry = 0.64 ± 0.15. Their data and the CCM extinction curve with 
the reported Ry value are shown in Figure 2. Motta et al. (2002) used a lens redshift of 
z = 0.83, found spectroscopically by Oscoz et al. (1997) using the Mg II absorption doublet 
and found Ry = 2.1 ± 0.9. If the extinction curve data given in Falco et al. (1999) are 
transformed to the same redshift as Motta et al. (2002) and normalized with the proper 
A(V), the data fall right on the curve suggested by Motta as seen in Figure 2. Falco et 
al. (1999) point out that using the technique of Jean & Surdej (1998) does not always lead 
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to accurate redshift and that using the wrong redshift will produce an incorrect extinction 
curve. Therefore, extinction curves produced from data with uncertain redshifts are not 
necessarily useful. 

Two important assumptions, made by all the studies listed in Table 1, are that one 
sightline is less reddened than the other, and that both sightlines have dust with the same 
extinction properties, i.e., the same value of Ry. If both images of the QSO are reddened 
but with different values of Ry then combining them to produce an extinction curve is not 
so straightforward. The wavelength dependence of dust extinction in any galaxy is unknown 
a priori, but for simplicity let's assume it follows the CCM relation. 

The CCM relationship at any wavelength, A, is: 

Let's assume that each QSO image has a dust column with different amounts and Ry values. 
Then, the magnitude difference between two lensed images is: 
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where Ry\ and ^4i(V) give the dust properties of the least attenuated, "unreddened" , image 
and Ry 2 and A 2 (V) of the most attenuated, "reddened", image (Wucknitz et al. 2003). 
Normalizing the extinction at any wavelength by the amount of extinction in the V-band 
gives: 
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Equation 7 gives the extinction curve that is actually measured for our simple example. It 
is the CCM relationship with an Ry value of 
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Therefore, the combination of two CCM extinction curves is also a CCM extinction curve 
but the resulting value of Ry is not an average of Ryi and Ry 2 . Thus, the measurement 
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of the difference between two lensed sightlines is actually the difference of two extinction 
curves, not an extinction curve itself. In order to avoid confusion, we propose the term "lens 
difference curve" be used for these measurements instead of labeling them extinction curves 
which can be misleading. 

As a quantitative example, in Figure 3 we show the lens difference curve for B0218+357, 
described by Munoz et al. (2004) with Ry = 12, fit to 6 photometric data points. We have 
reproduced this lens difference curve with a combination of different dust types for each 
image having Ry = 5.5 and 2.7. This result is not in any way unique. In fact, a very 
wide range of values of Ryu, can be produced by varying the values of Ry\ and Ry2- An 
example is shown in Figure 4, with A\{V) / A2{V) = 3. For this choice of Ai(V)/A2(V), we 
find 1.5 < Ry\2 < 56. Similar results are seen for other values of A 1 (V)/A 2 (V). It is even 
possible to produce negative values of Ry for such lens difference curves (Wucknitz et al. 
2003). 

In general, for any two QSO images, the value of A 1 (V)/A 2 (V) is unknown. This 
highlights the main difference between the pair method as applied to lensing galaxies and 
the pair method as applied to stars in the Galaxy and the Magellanic Clouds (Valencic et 
al. 2004; Gordon et al. 2003). In the latter case, the reddening, E(B-V), and the ratio of 
total-to-selective extinction, Ry, can be independently measured along the sightline to each 
star so that Ai(V)/A2(V) is a known quantity. Then, a unique extinction curve, within the 
measurement uncertainties, can be derived. Only in two limiting cases will the pair method 
produce an actual extinction curve for QSO images. First, if the wavelength dependence of 
the extinction is the same for dust along both sightlines then the resulting lens difference 
curve will be an extinction curve. In our simple CCM example, if the value of Ry is the 
same for the dust along both sightlines, then equations 7 and 9 produce a correct extinction 
curve and Ry value. Second, if A 2 (V) » ^-i(V) then it doesn't matter if the extinction 
parameters are different for the dust along the two sightlines. In this case, equations 7 and 
9 again produce the correct values when Ai(V) is set to zero. 

So the pair method for lensing galaxies can work but generally it is unknown which 
individual pairs of QSO images either have the same Ry, or have one image that is much more 
reddened than the other. Several hundred sightlines within the Milky Way have measured 
extinction curves (Valencic et al. 2004). The so-called average Milky Way extinction curve 
with Ry— 3.1 is, in fact, the typical extinction for diffuse sightlines in the Galaxy. It is not in 
any way an average curve for the Galaxy. Any two sightlines widely spaced in the Galaxy are 
likely to have very different values of Ry. The lensed QSO, B0218+357, where one sightline 
passes through a molecular cloud, is a case in point (Munoz et al. 2004). If the second 
sightline is lightly reddened then an accurate extinction curve for the dust in this molecular 
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cloud is being produced. But if the second sightline is also significantly reddened but with 
dust from the more diffuse interstellar medium in the intervening galaxy, then the resulting 
measurement is a lens difference curve not an extinction curve. This lens difference curve 
reflects the difference in the two sightlines' extinction curves and it is nontrivial to extract 
the original extinction curves from the lens difference curve. A very wide range of Ry values 
is listed for the lensing galaxies in Table f . Some fall in the normal range seen for Galactic 
sightlines, between 2.1 and 5.6, while other values are quite extreme. However, as can be 
seen from Figure 4, the normal Galactic values in Table 1 are no more likely to be correct 
than the extreme values. 

While it is not possible to extract the extinction curves from a lens difference curve, 
there is information about the dust in the lensing galaxy contained in the lens difference 
curve. While it may be possible to create a lens difference curve which follows the CCM 
relationship without the two individual extinction curves also following the CCM relation- 
ship, this seems unlikely. Therefore, if the lens difference curve can be fit with the CCM 
relationship regardless of the Ry value, this is evidence that the dust in the lensing galaxy 
follows the CCM relationship. In addition, if the 2175 A extinction bump is present in the 
lens difference curve, then there must be bump dust in the more reddened sightline at least. 
Munoz et al. (2004) focus on the extinction bump at 2175 A. Lens difference curves are 
shown for LBQS 1009-0252 and B0218+357, and the latter shows a bump at 2175 A which is 
reproduced here in Figure 3. Motta et al. (2002) also found the bump in their lens difference 
curve for SBS 0909+532. The presence of the bump in these lens difference curves implies 
that there is Milky Way- type bump dust, not just SMC- type dust along at least one of the 
component sightlines. 

Just because we do not see very large or very small values of Ry in the Local Group 
galaxies doesn't mean that exotic dust with very different properties can't exist in distant 
galaxies. Very simplistically, a small Ry implies that a greater mass fraction is contained 
in small dust grains and vice versa. However, extinction by large grains is less efficient 
than extinction by small grains because mass scales as a 3 while surface area scales as a 2 . 
Therefore, for any given gas abundance of grain materials and gas-to-dust-ratio, large grains 
will use a much higher fraction of the available material in order to produce the same amount 
of extinction. For instance models for dust along Galactic sightlines having high R V) such as 
the Orion sightline HD 37022 (Ry = 5.5), tend to run out of material to build large grains 
(Clayton et al. 2003). Therefore in general, sightlines with very small values of Ry are more 
likely to be physically possible than very large Ry values. 
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4. Conclusions 

• In general, the difference between two lensed sightlines in a galaxy does not produce 
an extinction curve, but a lens difference curve. It is not easy to decompose this lens 
difference curve into the two sightline extinction curves. 

• We show that the combination of two CCM type extinction curves having different 
values of Ry will produce a lens difference curve which follows the CCM relationship. 
The R v value of the lens difference curve is a function of the Ry values of the individual 
lensed sightlines and the ratio of their V band attenuations. 

• If one of the two QSO images is lightly reddened or if the dust along both sightlines 
has the same extinction wavelength dependence then the lens difference curve is an 
accurate extinction curve for the lensing galaxy. 

• Lens difference curves can only be accurately constructed for lenses with known redshift 
and magnification ratio, AM. If the redshift is wrong, the lens difference curve is 
wrong. 

• A small wavelength coverage results in lens difference curve CCM fits that are not well 
constrained. 

• The presence of a "bump" at 2175 A in a lens difference curve implies that there is 
Milky Way-like dust along at least one of the sightlines. 

• An accurate CCM fit to a lens difference curve, regardless of the fitted Ry values, 
implies that the two lensed sightlines have CCM-like dust extinction curves. 
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Table 1. Summary of Recent Work 
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a parenthesis indicate estimated redshifts with no spectroscopic confirmation 
b only available for objects with radio data 

"References.- (1) Falco et al. (1999); (2) Munoz et al. (2004); (3) Jean & Surdej 
(1998); (4) Motta et al. (2002); (5) Toft et al. (2000) 

d Motta et al. use spectroscopic data 

C AM determined from the best-fit extinction curve 
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Fig. 1.— Extinction data for B1600+434 (squares) fitted with a CCM curve with R v = 0.92 
(crosses) from Falco et al. (1999). Additional CCM extinction curves are plotted with 
Rvlow=0.65 and Rvhigh=1.3. This plot shows how weakly constrained the "fitted" CCM 
extinction curve is. The data are presented in rest wavelength space. 
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Fig. 2.— Motta et al.(2002) and Falco et al.(1999) data for SBS 0909+532 shown with CCM 
fits for R v = 2.1, the value obtained in the fit by Motta et al., and Ry = 0.64, the value 
reported by Falco et al. for their best-fit extinction curve. Also plotted are the Falco data 
using the Motta redshift and normalization, A(V). All data have been transformed into rest 
wavelength space using the reported redshifts. 
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Fig. 3.— Plot includes data from Munoz et al.(2004) for B0218+357, which fits a CCM 
(Rv = 12) extinction curve to photometric data in rest wavelength space. This plot shows 
two extinction curves with the following parameters: R V i = 5.5, R V 2 = 2.7, E(B—V) 1 = 0.7, 
E(B — V) 2 = 0.4 as well as the combined curve that would result which has Ryi2 ~ 12. 
Note that the choice of Ry\ and Ry2 are not unique. See text. 
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Fig. 4. — For any fixed ratio, A 1 (V)/A 2 (V), it is possible to obtain a wide range of combined 
Rv values, Ryi2- For this fixed value of A 1 (V) j 'A 2 (V) = 3, we find 1.5 < Rvn ^ 56. 



